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Abstract: A logarithm processing algorithm to measure beam transverse size and position is proposed and prelim¬ 
inary experimental results in Hefei Light Source II (HLS H) are given. The algorithm is based on only 4 successive 
channels of 16 anode channels of multianode photomultiplier tube (MAPMT) R5900U-00-L16 which has typical rise 
time of 0.6 ns and effective area of 0.8x16 mm for a single anode channel. In the paper, we hrstly elaborate the sim¬ 
ulation results of the algorithm with and without channel inconsistency. Then we calibrate the channel inconsistency 
and verify the algorithm using general current signal processor Libera Photon in low-speed scheme. Finally we get 
turn-by-turn beam size and position and calculate the vertical tune in high-speed scheme. The experimental results 
show that measured values ht well with simulation results after channel differences are calibrated and the fractional 
part of the tune in vertical direction is 0.3628 which is very close to the nominal value 0.3621. 
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1 Introduction 

In electron storage ring and collider complex, beam 
transverse profile measurement based on visible light 
from synchrotron radiation (SR) is widely used [l|. In 
general, CCD can be applied to measure beam size and 
position steadily and online. But conventional CCD has 
a long exposure time and very low frame rate. In HLS 
II, there has been installed a CCD camera to measure 
beam profile and monitor the machine status. The CCD 
camera has the maximum frame rate of 205 Hz at the full 
resolution (640x480), and the electron beam go through 
about 22k turns within single frame interval. As a re¬ 
sult, the measurement results are indeed the average of 
multiple bunches in multiple turns. 

In order to track the beam size and position of single 
bunch, it is necessary to acquire transient light intensity 
information utilizing high speed photodetector and data 
acquisition system to ensure that the transformed light 
pulse are not piled up along the pulse train. The typical 
schemes are VEPP-4M Q and CESR-c ^ which calcu¬ 
late the turn-by-turn or bunch-by-bunch beam size and 
position to study beam-beam interaction, beam dynam¬ 


ics, etc. In these measurement schemes, dozens of chan¬ 
nel information up to 16 or 32 are sampled in parallel to 
fit the shape of Gaussian beam spot, which increase the 
costs and complexity. Therefore, we put forward a log¬ 
arithm processing algorithm to calculate the beam size 
and position which needs only 4 successive channels, and 
carry out some experiments to verify the algorithm. 

2 Principle of the measurement system 

At the moment there are 4 applications in the visi¬ 
ble diagnostic beam line at 15° outlet of bending mag¬ 
net in HLS 11. Eig. [T] shows the block diagram of the 
measurement system for CCD and vertical MAPMT. 
Synchrotron radiation is first reflected by water-cooled 
Oxygen Eree Copper (OEC) mirror in vacuum chamber 
and then the visible light pass through several mirrors 
and a main focusing lens with magnification of I be¬ 
fore reaching the optical table. MAPMT and CCD Q 
share the same magnification of 2 by front-end optical 
path. Then the beam spot reaches the horizontal and 
vertical MAPMT cathode with magnification of 6 and 
magnification of 24 respectively. Magnification can be 
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altered to meet the need in different machine status. 
The two MAPMTs are fixed on three-dimensional trans¬ 
lation stage so as to find the image plane conveniently. 
Front-end circular slit constraints the opening angle of 
SR within ±3.8 mrad to decrease the depth of field error 


and diffraction error @ . A band pass filter with 500 nm 
peak wavelength and 10 nm full width at half maximum 
(FWHM) is used to allow quasi-monochromatic light to 
get to the sensitive area of MAPMT. 


source point 



The measurement system works in two schemes. In 
low-speed scheme, photocurrent of 4 successive chan¬ 
nels is sampled by general current signal processor Lib¬ 
era Photon from Instrumentation Technologies and then 
sent to remote PC by EPICS/C A protocol, where we 
can calibrate the channel inconsistency, implement the 
logarithm processing algorithm to extract beam size and 
position, and control the negative high voltage to pro¬ 
tect the MAPMT. In high-speed scheme, photocurrent is 
firstly conditioned by analog front-end electronics (AFE) 
including current to voltage conversion, preamplifier and 
low pass filter. Then transformed signals are sampled si¬ 
multaneously by high speed oscilloscope DSO9104A and 
acquired over Local Area Network (LAN) by remote PC 
where we can calculate turn-by-turn beam size and po¬ 
sition, perform East Eourier Transformation (EET), etc. 

3 Logarithm processing algorithm ^ 


in Eig.[2l The rectangular region of A, A, A and A rep¬ 
resents effective photosensitive areas with 0.2 mm gap 
between them. 



Photocurrent of A, A, A and A, which derives from 
4 successive channels, is proportional to the integral of 
light intensity and can be expressed by the subtraction 
of two error functions shown in Eq. ©• 


As SR light irradiated by electron beam from bend¬ 
ing magnet in storage ring presents Gaussian distribu¬ 
tion, we assume SR light intensity distribution function 
4>(x) as Eq. ([T]). 


4>(x) = 4>oexp 



(1) 


where 4>o is the maximum light intensity; S is beam dis¬ 
placement that center of the beam spot offsets specified 
center of the MAPMT; a is beam size. 

Based on the assumption of Gaussian distribution, 
the SR light intensity distribution at MAPMT is shown 


A = 
A = 
A = 
A = 


XI 

^Xl 


-X2 

— X3 


^(x)dx = A 
4>(x)dx = A 
^(x)dx = A 
^(x)dx = A 


erf -erf (^) 


x/2 


^/2 


erf(Z^)_erf(Z^) 
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^/2 


erf( ±'^ )-erf( ± 
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where erf(x) = ^ f^e~^^dt; Xi=0.1 mm, ^2=0.9 mm, 
X 3 =l.l mm and ^4=1.9 mm; A is the maximum anode 
photocurrent. 
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3.1 Ideal case 


In ideal condition when the 4 continuous channels 
have the same response characteristic, we define ideal 
size signal 5'in(cr,(5) and ideal position signal Pin(cr,(5) in 
Eq. (|3]). 


' S,„(a,6) = [\n{l2ls/Iil4)r''\ 


( 3 ) 


Substituting the integral result shown in Eq. m into 
Eq. (|3]) and simplifying it, we can obtain how the size and 


position signal vary from size and position illustrated in 
Fig. El When beam size increases from 0.2 mm to 2 mm 
with 0 mm beam displacement, size signal 5'in(cr,^^) has 
a good linear relation with size a and the linear fitting 
equation can be expressed by Sin{<j) = 0.03924 + 0.69cr. 
When beam size varies within 0.8mm ^ 2mm, beam 
displacement S has an effect on normalized size signal 
within 1%. Meanwhile, position signal Pi^{a,S) has a 
good linear relation with beam displacement when beam 
displacement varies from -2 mm to 2 mm and linear fit¬ 
ting equation is P\n{S) = 2S. Eurthermore, beam size has 
an impact on position signal sensitivity within 1%. 




Fig. 3. The simulation results of ideal beam size and position signal in different conditions: (a) size signal versus 
beam size when (5=0 mm; (b) normalized size signal versus beam displacement when a = 0.8 ~ 2.0mm; (c) position 
signal versus beam displacement when cr = 0.82.0mm; (d) position signal sensitivity versus beam size. 


Erom the above simulation results, we can get the 
approximate equations expressed in Eq. 

f (CT, (5)« S'ln (a) = 0.03924 + 0.69a, 
1Pi„(a,5)«P,„(a) = 2(5. ^ ^ 


3.2 Non-ideal case 


Given the channel inconsistency resulting from anode 
uniformity and cross-talk among channels, we introduce 
the channel gain factor = 1, 2, 3, 4) to correct chan¬ 
nel inconsistency based on linear approximation, as is 
shown in Eq. ®. 


5,n(a,<5) 

an(cT,(5) 


ln(PJ3/P/4) 


- 1/2 


[ln(fl253M<74)+ln(/2/3//l/4)]“'^t 


Mhh/hh) 

ln(glfl2/fl3g4)+ln(JlJ2/l3l4) 
Ill((72fl3/fll54) +lll(/2l3/Pl4) 


( 5 ) 


where Pin(cr,(5) is actual size signal, Pin(cr,(5) is actual 
position signal, = gili{i = 1, 2, 3, 4) represents ac¬ 
tual anode photocurrent from 4 successive channels with 
channel inconsistency. 

If we define d, = In(g2fl3/5ifl4), dp = ln(gi52/5354) and 
put Eq. (E]) into Eq. (EJ, we can get: 


' 5,n(a,,5)«5,„(a) 

l + J,S'f„(a) 

< Pln(a,(5)RiPln((5) 

l + 4S'f„(a) 



l + dsSf^{a 


where dg and dp are gain correction factors. 

According to Eq. (|6]), we know that actual position 
signal has a good linear relation with beam position, but 
actual size signal presents a huge nonlinearity with size. 
In addition, dg and dp can be induced by slope and in¬ 
tercept of Pin(cr,(5) if beam size and magnification can be 
acquired in advance. Hence, modified size signal Pin(cr,(5) 
and modified position signal Pin(cr,(5) can be expressed in 
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Eq. (j7|). 


su^,s) 

< 


-An(cr,(5) 


ln{iJs/hh)-l 

Inihh/W-d^ 



\n{his/ij4)-l' 


( 7 ) 


4 Experimental results 

4.1 Low-speed scheme 

In low-speed scheme, the horizontal measurement re¬ 
sults are described. In order to simulate the movement 
of beam center, we move the three-dimensional trans¬ 
lation stage for horizontal MAPMT horizontally. The 
movement range is from 13.0 mm to 16.8 mm with 0.1 
mm step. The process of movement is very short, so the 
beam size measured by CCD remains unchanged. 

Figure |4] shows how the average anode output cur¬ 
rent varies with the moved position. CH4, CH5, CH6 and 
CH7 are corresponding to channel number 4, 5, 6 and 7 
respectively. When the horizontal translation stage in¬ 
creases, average anode current of CH4 and CH5 increases 
and it is just opposite for CH6 and CH7, which roughly 
indicates the beam spot center is moving from CH7 to 
CH4. In the middle point of 14.9 mm, anode current of 
CH4 and CH7 are approximate the same, and so does it 
for CH5 and CH6. Thus this position can be thought to 
be the absolute center. 



Fig. 4. Average anode current of 4 successive channels. 

Putting average anode current into Eq. (j5j), we can 
get how the actual position signal varies with the simu¬ 
lated beam displacement shown in Fig. O The horizon¬ 
tal coordinate axis has been adjusted by subtracting the 
absolute center. Applying linear fitting to the raw dis¬ 
crete points, we can get the fitting equation with slope of 
1.7894 and intercept of-0.1206. Actual position signal 
has a linear relation with simulated beam displacement 
as described in Eq. However, the fitting curve does 
not cross the point (0, 0), and slope does not equal to the 
simulation result in ideal case shown in Fig. ^c). The 
difference mainly comes from the channel inconsistency 
and must be calibrated. 



Fig. 5. Actual position signal with channel inconsistency. 

Substituting the magnifying beam size imaged at the 
horizontal MAPMT cathode into Eq. ([6j), we can get the 
measured gain correction factors =0.0590 and dp=- 
0.0675 on the basis of the fitting equation shown in Fig.[5l 
As a consequence, the relationship between the modi¬ 
fied position and simulated position can be illustrated 
in Fig. [6l Linear fitting results show that the modi¬ 
fied position signal has a good linear relation with simu¬ 
lated beam displacement and the fitting curve cross point 
(0, 0) with slope of 2.002 and intercept of -0.008 which 
agrees well with the simulation result in Fig .3(c). 



Fig. 6. Modified position signal. 

Beam size calculated by direct equation in Eq. ^ 
and modified equation in Eq. 0 is illustrated in Fig. [3 
respectively. Direct and modified beam size gradually 
decreases when beam center is moved from -1.8 mm to 
1.9 mm although the modified beam size signal which 
is proportional to beam size should remain almost un¬ 
changed according to Fig .3(b). But the modified beam 
size is more close to the nominal value 0.332 mm mea¬ 
sured by CCD and is more reliable. 



Fig. 7. Direct and modified beam size. 


4.2 High-speed scheme 

In high-speed scheme, we take the vertical measure¬ 
ment results for example. Single bunch fills in the stor¬ 
age ring with 3 mA beam current and electron beam is 


010201-4 





















Submitted to “Chinese Physics C' 


excited by swept frequency signal which comes from the 
tracking generator module in the spectrum analyzer. Sig¬ 
nals from AFE are continuously sampled 500,000 points 
with 10 GS/s sample rate for 4 channels. External trigger 
signal with frequency of 4.534 MHz synchronized with 
the revolution frequency is applied to keep pace with the 
pulse signal produced by SR light. Eig.[8]shows the peak 
values of each turn extracted from raw data points with 
peak search algorithm. As the sample points per one 
period is about 2200, we can sample 227 turns data in 
total. CHS, CH9, OHIO and CHll are corresponding to 
channel number 8, 9, 10 and 11 respectively. Peak value 
of CH9 and OHIO is larger than CHS and CHll indi¬ 
cating that beam center is located between the CH9 and 
CHIO which satisfies the conditions of the algorithm. 

-0.15 
-0.20 
-0.25 

§ -0.30 

O) 

S -0.35 

o 

^ -0.40 

03 

Q. -0.45 
-0.50 
-0.55 

0 40 80 120 160 200 240 

turn number 

Fig. 8. Peak values of each turn for 4 successive channels. 

Applying FFT to all the 500,000 points raw data 
and enlarging the view of spectrum of CH9 in frequency 
range of 0 Hz to 10 MHz shown in Fig. [9l we can clearly 
see vertical synchronous oscillation sidebands appear at 
frequency of 1.620 MHz, 2.900 MHz, 6.160 MHz and 
7.440 MHz except for two main harmonic components 
at the frequency of 4.530 MHz and 9.060 MHz which 
are roughly equal to the integral multiple of revolution 
frequency. These sideband elements result from the in¬ 
teractions between the synchronous light and slit. In 
HLS H, frequency at 1.620 MHz and 6.160 MHz cor¬ 
respond to betatron oscillation in vertical direction, and 
frequency at 2.900 MHz and 7.440 MHz are the mirrored 
components of vertical betatron oscillation. 



Fig. 9. Spectrum of CH9. 


With the channel inconsistency being calibrated and 
logarithm processing algorithm being adopted with 227 
turns peak values, we can get the turn-by-turn beam size 
and position shown in Fig. [TOl From Fig. flOl a). beam 


size varies between 64.5 /im and 93.5 /im with the aver¬ 
age of 76.4 /im. From Fig. flOlb). beam position varies 
between -8.2 /im and 22.6 /im and the average value is 
8.0 /im . Meanwhile, we record the beam size with CCD 
which tells that the beam size is about 70 jam. Our aver¬ 
age result of beam size is larger than the beam size from 
CCD for about 8% which may result from the amplifica¬ 
tion of light path. 



0 40 80 120 160 200 240 

turn number 

Fig. 10. Beam size and position for 227 turns: (a) 
beam size; (b) beam position. 

Applying FFT to the turn-by-turn beam position and 
getting rid of the dc component, we can obviously see a 
peak signal at the frequency of 1.645 MHz correspond¬ 
ing to the vertical betatron oscillation sideband shown 
in Fig. [TTJ So the vertical fractional tune is 0.3628 with 
4.534 MHz revolution frequency. The measured result is 
in good agreement with the nominal value 0.3621. 



Fig. 11. Spectrum of turn-by-turn position. 


5 Conclusion 

A logarithm processing algorithm used to calculate 
beam transverse size and position based on SR light in 
HLS II is developed and preliminary experimental results 
are given. The measurement results are in good agree¬ 
ment with the simulation results after the channel incon¬ 
sistency is calibrated. Vertical fractional tune of 0.3628 
is obtained according to the turn-by-turn beam position 
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which is very close to the nominal value 0.3621. More 
experiments will be conducted in multiple bunches filling 
mode, and effective ways to excite the obvious change of 
beam size should be applied to study beam instability 
and beam dynamics. 
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